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Abstract. High-resolution non- ideal magnetohydrodynamical simulations of the turbulent mag- 
' netized ISM, powered by supernovae types la and II at Galactic rate, including self-gravity and 
non-equilibriuim ionization (NEI) , taking into account the time evolution of the ionization struc- 
ture of H, He, C, N, O, Ne, Mg, Si, S and Fe, were carried out. These runs cover a wide range 
(from kpc to sub-parsec) of scales, providing resolution independent information on the injection 
scale, extended self-similarity and the fractal dmension of the most dissipative structures. 
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r 1. Introduction 

O ,. In star forming disk galaxies, matter circulation between stars and the interstellar 
medium (ISM), in particular the energy input by supernovae, determines the dynamical 
and chemical evolution of the ISM, and hence of the galaxy as a whole. So far ISM models 
used radiative cooling assuming collisional ionization equilibrium (CIE), which is a good 
approximation providing the cooling time is much longer than the recombination times 
(e.g., Kafatos 1973). A condition verified for most of the ions for temperatures > 10 5 ' 8 
— . K. At lower temperatures departures from equilibrium are expected. While in CIE the 
ionization fractions depend only on the temperature and are sharply peaked, in NEI 
these same fractions depend on the dynamical and thermal history of the plasma. These 
departures affect the local cooling, which is a time-dependent process that controls the 

. flow dynamics, feeding back to the thermal evolution by a change in the density and 
internal energy distribution, which in turn modifies the thermodynamic path of non- 

. equilibrium cooling (Breitschwerdt & Schmutzler 1999). 



2. ISM Modelling 

Following Avillez & Breitschwerdt (2007) we study the evolution of the ISM in a 
patch of the Galaxy with (x,y) ^ 1 kpc size in the Galactic plane, and \z\ ^ 10 
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kpc perpendicular to it using non-ideal MHD equations coupled to the time-dependent 
calculation of the ionization structure of the plasma. The adopted model includes: (i) 
SNe types la and II occurring at the Galactic SN rate - 40% of the SNe II occur in the 
field (ii) gravitational field of the stellar disk, (iii) local self-gravity, (iv) heat conduction, 
(v) backgroung heating due to the UV photon field, (vi) mean (of 3.0/1 G) and turbulent 
magnetic field components corresponding to a total field of 4.5 [i G, and (vii) radiative 
cooling calculated on the spot and using solar abundances of Asplund et al. (2005). The 
NEI calculation takes into account the time-dependent ionization structure of the 10 
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most important elements in nature - H, He, C, N, O, Ne, Mg, Si, S, and Fe - and includes 
collisional ionization by thermal electrons, autoionization, charge exchange reactions, 
radiative and dielectronic recombination, photoionization, comptonisation, Auger effect 
and ionization of H and He by photo- and Auger electrons. 



3. Fractal Dimensions 




Figure 1. Left: The scalings C(p)/C(3) vs. order p. Black, red and blue lines refer to Kolmogorov 
(1941), She-Leveque (1994) and Boldyrev (2002) models, respectively; Bullets refer to data of 
Benzi et al. (1993); Triangles (black and red) and green squares refer to CIE HD (1.25 pc and 
0.625 pc resolutions) and MHD (1.25 pc resolution) runs of Avillez & Breitschwerdt (2007); Blue 
circles correspond to NEI MHD (0.5 pc resolution) run discussed here. Right: The Hausdorff 
dimension D of the most dissipative structures vs. C(p)/C(3), with p — 2, 4, 10 shown in the 
left panel. 

We measured the velocity structure functions (Svf) and determined the correspond- 
ing scalings C(p)/C(3) and the Hausdorff dimension D of the most dissipative structures 
(Figure [1]). D oscillates around 2, indicating that turbulent energy is dissipated preferen- 
tially through shocks in the HD cases and a combination of shocks and current sheets in 
the MHD runs. In the presence of the magnetic field (ideal and non- ideal calculations) D 
departures from 2 due to the anisotropy induced by the magnetic field. The NEI and CIE 
runs have similar fractal dimensions for the most dissipative structures and reproduce 
experimental and numerical results previously published. 
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